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Charge spin separation, pseudogap formation and phase diagrams are studied in two and four 
site Hubbard clusters using analytical diagonalization and grand canonical ensemble method in 
a multidimensional parameter space of temperature, magnetic field, on-site Coulomb interaction 
(U > 0), and chemical potential. The numerically evaluated, exact expressions for charge and 
spin susceptibilities provide clear evidence for the existence of true gaps in the ground state and 
pseudogaps in a limited range of temperature. In particular, Mott-Hubbard type charge crossover, 
spin pseudogap and magnetic correlations with antiferromagnetic (spin) pseudogap structure for 
two and four site clusters closely resemble the pseudogap phenomena and the normal-state phase 
diagram in high T c superconductors. 
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Understanding the effects of electron correlations and 
pseudogap phenomena Q, 0, 0, 0, H, El in the cuprate 
superconductors comprising of many different phases is 
regarded as one of the most challenging problems in con- 
densed matter 0. Although the experimental deter- 
mination of various inhomogeneous phases in cuprates 
is still somewhat controversial |8j, the underdoped high 
T c superconductors are often characterized by crossover 
temperatures below which excitation pseudogaps in the 
normal-state are seen to develop [9(- There is also com- 
pelling evidence for the existence of quantum critical 
points (QCPs) in underdoped 0, 0, and opti- 
mally doped materials ^3 as observed in resistivity mea- 
surements in Nd2- a; Ce :! ;Cu04± ( 5, P^-^Ce^CuO^-a and 
La2- a; Sr a ;Cu04. 

The charge-spin separation 0, clearly identifi- 
able Mott-Hubbard (MH), antiferromagnetic and spin 
crossovers contain generic features which appear to be 
common for small clusters and large thermodynamic sys- 
tems 16]. Studies of conventional Mott-Hubbard and 
magnetic phase transitions 0, 0, 0, |2?J |2l| have 
tended to concentrate on macroscopic systems contain- 
ing a large, and essentially an infinite number of parti- 
cles. In attempts to address some of the above, the Hub- 
bard model has been discussed within the exact Lieb-Wu 
(LW) equations in one dimension (Id) [22l I23I |24| and 
a wide variety of approximation schemes in higher di- 
mensions [2^, |2(| I23. Most theories originating from 
the Bethe-anstaz, such as LW, involve coupled nonlin- 
ear integral equations that have to be solved numerically 



for every set of parameters. Numerical uncertainties as- 
sociated with such solutions severely limit their appli- 
cations when calculating subtle features at intermediate 
values of temperature and other parameters. Although 
some properties of Hubbard clusters have been calcu- 
lated Hi IH |U IH HI |H H3, many questions re- 
main with regard to microscopic origins of charge-spin 
separation and pseudogap behavior, short range correla- 
tions and weak singularities (crossovers) at finite tem- 
perature |35ll36|. 

In this work, we study the phase diagrams for the two 
and four site Hubbard clusters 0] using analytical diag- 
onalization combined with the grand canonical ensemble 
in a multidimensional parameter space of temperature 
T, magnetic field h, on-site Coulomb interaction U > 0, 
and chemical potential [i (with hopping parameter t = 1). 
Our calculations for finite clusters are based on the ex- 
act analytical expressions for the eigenvalue E n of the 
n th many-body eigenstate, grand partition function Z 
(where the number of particles N and the projection of 
spin s z are allowed to fluctuate) and its derivatives with- 
out taking the thermodynamic limit and hence involve 
no approximations. The grand canonical potential Qu 
for interacting electrons is 



-Tin £ , 

n<N H 



(1) 



where N n and s z n are the number of particles and the 
projection of spin in the n th state respectively. The di- 
mension Nh of the Hilbert space in depends on the 
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FIG. 1: Temperature T vs chemical potential n phase dia- 
gram for the two site cluster at U = 6, h = 0. The antiferro- 
magnetic pseudogap phase is shaded and bounded above by 
Taf(h)- It vanishes in regions I and II, which are charge bi- 
furcation phases representing (pseudo) charge gaps at T > 0. 
The charge and spin crossover temperatures T c (n) and T a (fi) 
are obtained from maxima in charge and spin susceptibilities. 
Since there is particle-hole symmetry, only [i values above 
half-filling, i.e. n > /j,q, are shown. Note that in regions I and 
II, there is strong charge-spin separation. 

number of sites, satisfying Njj = 4 2 for the two site clus- 
ter, and Nh = 4 4 for the four site cluster. The electron 
charge susceptibility Xc(m) or the corresponding thermo- 
dynamic density of states, p(fi) = , describes the 
local spectral characteristics of charge excitations and 
fluctuations in the number of electrons 

(N 2 )-(N) 2 = T^^. (2) 

O/i 

The spin susceptibility Xs(m) or s P m density of states 
cr(jLi) = dt ^ h ^ describes the local spin excitations as pa- 
rameters ft, /i or T are varied, where spin fluctuations 
{'(As 2 ) 2 ) closely follow the variation of spin susceptibil- 
ity x with respect to /i and h 

((^) 2 )-<* 2 > 2 =T^p. (3) 

Once all the many-body eigenvalues of the Hubbard 
clusters are known, it is straightforward to calculate the 
ground state properties and these results are reported 
elsewhere [lfj. Using the same eigenvalues, we have eval- 
uated the exact grand partition function and thermal av- 
erages such as magnetization and susceptibilities numer- 
ically as a function of the set of parameters {T, h, [i, U}. 



FIG. 2: Temperature T vs chemical potential /j, phase diagram 
for the four site cluster at U = 6, h = 0. Regions I and II 
are quite similar to the ones found in the two site cluster, 
again showing strong charge-spin separation. However, note 
the (new) mixed phase which consists of charge, spin and 
AF pseudogaps, that is present in this cluster. The inset 
shows the same diagram with a larger scale for the crossover 
temperature. Labels of crossover temperatures are suppressed 
for clarity. 



Using maxima and minima in spin and charge suscepti- 
bilities, phase diagrams in a T vs fi plane for any U and 
h can be constructed. 

Among many interesting results rich in variety for 
U > 0, sharp transitions are found between phases with 
true charge and spin gaps in the ground state; for in- 
finitesimal T > 0, these gaps are transformed into 'pseu- 
dogaps' with some nonzero weight between peaks (or 
maxima) in susceptibilities monitored as a function of 
doping (i.e. /i) as well as h. We have also verified the 
well known fact that the low temperature behavior in the 
vicinity of half-filling, with charge and spin pseudogap 
phases coexisting, represents an antiferromagnetic insu- 
lator in the Hubbard clusters |16(. However, away from 
half Riling, we find very intriguing behavior in thermo- 
dynamical charge and spin degrees of freedom. At low 
temperature, new peak structures in the charge Xc(m) 
and zero magnetic field spin Xs{[*) susceptibilities are 
observed to develop between two consecutive peaks, 
there exists a pseudogap in charge or spin degrees. Open- 
ing of such distinct and separated pseudogap regions at 
low temperature, the signatures of corresponding charge 
and spin separation, is seen in both the two and four site 
clusters, away from half filling. 

In Fig. ^ we show the phase diagram for the two site 
cluster at U = 6 and h = 0, where T c (fi) describes the 
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high temperature, MH insulator-metal transition as fi 
varies. At low temperature, this dependence bifurcates 
and T c i(ji) and T C 2(/x) denote the crossover temperatures 
at lower and higher doping (respectively) corresponding 
to similar MH like transitions. The low temperature 
regions I and II, bounded by these crossover tempera- 
tures, define a charge bifurcation phase which persists 
up to a high doping level, [i c i. In this charge pesudo- 
gap, there are strong spin fluctuations centered around 
fi s at low temperature which decay when either fi — > [i c \ 
or /i — > /i C 2 (i.e. near the boundaries where charge gaps 
melt). In the same figure, the antiferromagnetic crossover 
temperature Taf(^) is obtained by monitoring the spin 
susceptibility peaks as a function of the applied field at 
a given chemical potential /i, as seen in measurements 
of Neel temperature doping dependence 0. If h c de- 
notes the field value corresponding to the peak closest 
to zero, Taf is defined as the temperature at which the 
critical field h c (/j) for the onset of magnetization van- 
ishes (melting the AF spin gap). Taf can also obtained 
from zero field staggered spin susceptibility. The spin 
crossover temperature T a (fi), associated with the open- 
ing of the zero-magnetic-field spin pseudogap in Fig. ^ 
denotes the temperature below which strong spin pseudo- 
gap correlations are observed to develop. This has been 
suggested as a precursor to superconductivity [37| ■ Note 
that spin fluctuations are significantly suppressed below 
T s and how these crossover temperatures and critical /x 
values define various regions such as, a) charge pseudo- 
gap, b) AF pseudogap, c) spin pseudogap and d) normal 
phases. 

The four site phase diagram, Fig. [21 with U — 6 and 
h = 0, appears to be more complex with several charge 
bifurcation phases above (or below) half filling, start- 
ing at the electron (hole) doping level near 1/8 filling. 
However, there is a clear self-similarity in the right-most 
bifurcation region for the 2-site and the 4-site clusters 
(Figs, ^and [2J. This bifurcation region retains strong 
charge pseudogap stability as in the two site cluster. We 
believe that these regions of phase space, with strong 
spin and charge fluctuations, are quite relevant to the 
high T c cuprates and other similar materials, where dop- 
ing of electrons or holes introduces dramatic changes in 
their physical properties. 

We have followed the behavior of pseudogaps in this 
region as a function of the on-site Coulomb repulsion U. 
With increasing U, several features are observed to de- 
velop; (a) separation of charge and spin boundaries away 
from half-filling, (b) opening of a pseudo charge gap, 
(c) large spin fluctuations inside this charge gap region. 
This phase with a charge gap closely resembles the inho- 
mogeneous commensurate phase found in Ref.0, where 
the mobile holes in the quasi-one dimensional structures 
(stripes) acquire a spin gap in spatially confined Mott- 
insulating regions @, ■ 

The doping dependence of the normal-state spin pseu- 
dogap shows that at low temperature, it is stable in the 
underdoped regime, persists at optimal doping and dis- 



appears in the overdoped regime at /i = fj, s . In the phase 
diagrams in Figs. ^ and [21 at zero temperature, we no- 
tice several quantum phases and corresponding quantum 
critical points (QCPs), at fj, = fj, c i, [x s , /x C 2. A compari- 
son of the two phase diagrams for the two and four site 
clusters, in Figs. ^ and [21 reveals many common features 
in addition to the QCPs. As /i increases, there is a sharp 
transition in both clusters at fi c i from a Mott-Hubbard 
antiferromagnetic insulator into a phase with charge and 
spin separation and gaps. A similar behavior has been 
observed experimentally [THIT^] in the cuprates. At crit- 
ical doping /i c i, the antiferromagnetic phase disappears 
and at higher doping, in regions I and II, the spin and 
charge pseudogap phases coexist with one another inde- 
pendent of how strong U is. These two regions are sepa- 
rated by a boundary where the spin gap vanishes. At crit- 
ical doping /x s with T s — > 0, the zero spin susceptibility 
Xs at zero temperature exhibits a sharp maximum. Thus 
the behavior of the critical temperature T s (n), which falls 
abruptly to zero at critical doping fi s , implies |8Lll0j| that 
the pseudogap can exist independently of possible super- 
conducting pairing. As mentioned above, as T — > and 
H = ji Sl the spin gap disappears while the charge gap 
prevails up to /i = [i c i . Up on further doping, the charge 
gap vanishes at fi C 2 and Fermi liquid behavior is restored 
due to full charge-spin reconciliation. Figs. ^ and [21 are 
consistent with the existence of pseudogap phases and 
quantum phase transitions at fj, s in the high T c super- 
conductors, when the ground state spin gap disappears 
H H El The QCPs separate the spin pseudogap 
phase from the quantum spin liquid state, coexisting with 
the charge pseudogap. We have also seen that a reason- 
ably strong magnetic field has a dramatic effect (mainly) 
on the QCP at fx s . 

The spin susceptibility Xs i n the underdoped pseu- 
dogap region I, defined by [i c \ < [l < jjL s , displays an 
anomaly at temperature T s below which the spin de- 
grees of freedom become suppressed. The gapped spin 
excitations imply that spin singlet states which exist be- 
tween T s and T e \ in Fig. ^ ca n be considered as pre- 
formed pairs with properties different from FL. As tem- 
perature decreases, approaching T c \ from above, Xs de- 
creases while Xc increases due to strong charge fluctua- 
tions in the vicinity of T c ± signaled by a sharp peak in 
the excitation spectrum, consistent with ARPES mea- 
surements 39]. The thermal fluctuations close to the 
boundaries T cl of the charge bifurcation phase destroy 
the charge pseudogap and provide conditions at which 
the spin pseudogap state coexists with a charge liquid 
background. 

In Figs, ^and [21 we also show the crossover tempera- 
ture Taf(ij-) deduced from maxima in the antiferromag- 
netic spin fluctuations cr(fi). The phase diagram for the 
two site cluster shows an island of stability for the antifer- 
romagnetic phase and the corresponding crossover tem- 
perature, Taf{^), decreases monotonically as the chemi- 
cal potential fj, increases. In the underdoped regime close 
to half filling, the antiferromagnetic phase in the two site 
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2-site phase diagram: U=0 



4-site phase diagram: L'=0 
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FIG. 3: The single particle or 'noninteracting' (U = 0) case, 
illustrating the charge and spin, peaks and valleys for the two 
and four site clusters. Note how the charge and spin maxima 
and minima in dos follow one another indicating the absence 
of charge-spin separation. 



cluster is fully separated from the charge bifurcation re- 
gion. In contrast, the behavior of Taf with increasing 
fi is nonmonotonic for the four site cluster and there is 
some overlap (near Taf ~ T c \) in the underdoped regime 
between the antiferromagnetic phase and a charge bifur- 
cation region. At higher doping within the (right-most) 
bifurcation regime (i.e. regions I and II), Taf = and 
the charge-spin separation is very similar to that seen in 
the two site cluster. 

As an important footnote, we show how the charge 
and spin peaks (and valleys) follow one another for both 
the two and four site clusters when U = in Fig. [3] 
(in sharp contrast to the U = 6 cases, in regions I and 
II where charge (as well as spin) maxima and minima 
are well separated) indicating that there is no charge- 
spin separation here. In this limiting, single particle case 
for the two site cluster, there are no temperature driven 



charge bifurcation phases away from half filling and hence 
the corresponding T cl = 0. In the entire range — 2t < 
fi < 2t, the charge and spin fluctuations directly follow 
one another, i.e. T c \ = T C 2 = T s without charge-spin 
separation. 

In summary, we have illustrated the phase diagram 
and the presence of temperature driven crossovers, quan- 
tum phase transitions and charge-spin separation for any 
U 7^ in the two and four site Hubbard clusters as 
doping (or chemical potential) is varied. Our bottom- 
up approach and exact thermodynamics for small clus- 
ters, when monitored as a function of doping, displays 
the presence of clearly identifiable, temperature driven 
crossovers into new phases and distinct transitions at 
corresponding QCPs in the ground state, seen in large 
thermodynamic systems. It appears that the short-range 
correlations alone are sufficient for pseudogaps to form 
in small and large clusters, which can be linked to the 
generic features of phase diagrams in temperature and 
doping effects seen in the high T c cuprates. The pseu- 
dogap features and the variations of \c Xs, Taf with /i 
as well as the existence of QCPs suggest that the nor- 
mal state spin singlet pseudogap, closely linked to short 
range correlations, can also exist in small clusters. The 
exact cluster solution shows how charge and spin gaps 
are formed at the microscopic level and their behavior 
as a function of doping (i.e. chemical potential), mag- 
netic field and temperature. The pseudogap formation 
can also be associated with the condensation of spin 
and charge degrees of freedom below (spin and charge) 
crossover temperatures T s and T c \ (respectively). Fi- 
nally, the two and four sites clusters share very important 
intrinsic characteristics with the high T c superconductors 
apparently because in all these 'bad' metallic high T c ma- 
terials, local interactions play a key role. 

This research was supported in part by the U.S. De- 
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